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Within first-order phase transitions, we investigate the pre-transitional effects due to the nonper- 
turbative, large- amplitude thermal fluctuations which can promote phase mixing before the critical 
temperature is reached from above. In contrast with the cosmological quark-hadron transition, 
we find that the rapid cooling typical of the RHIC and LHC experiments and the fact that the 
quark-gluon plasma is chemically unsaturated suppress the role of non-perturbative effects at cur- 
rent collider energies. Significant supercooling is possible in a (nearly) homogeneous state of quark 
gluon plasma. 
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I. INTRODUCTION 

It is possible to model the gross general features of a 
phase transition from a quark-gluon plasma (QGP) to 
a hadronic phase through a phenomenological potential 
with a scalar order parameter [jjj. Assuming the tran- 
sition to be discontinuous, or first-order, as suggested 
by some recent lattice QCD simulations g, the QGP is 
cooled to a temperature T± , where a second minimum ap- 
pears, indicating the presence of a hadronic phase. With 
further cooling, the two phases become degenerate at the 
critical temperature T c , with a free energy barrier which 
depends on physical parameters characterizing the sys- 
tem, such as the surface tension (a) and the correlation 
length (£). This general behavior models both the cos- 
mological quark-hadron phase transition and the produc- 
tion of a QGP during heavy-ion collision experiments, as 
those under way at the RHIC and planned for the LHC. 
In the latter case, the plasma generated by the collision 
expands and cools, relaxing back to the hadronic phase. 
Recent interest has been sparked by the possibility that 
this relaxation process is characterized by the formation 
of disoriented chiral condensates (DCC), which are co- 
herent pion condensates similar to the domains typical 
of quenched ferromagnetic phase transitions The 
nonequilibrium properties of this relaxation process and 
DCC formation has also been studied as a first order 
chiral phase transition where the supercooled phase may 
naturally lead to a "quenched" initial condition ||. 

Recent work on the dynamics of weak first-order phase 
transitions have shown that, in certain cases, it is possi- 
ble to have nonperturbative, large-amplitude fluctuations 
before the critical temperature is reached, which promote 
phase mixing j(| . Studies performed in the context of the 
cosmological electroweak phase transition |24j and quark- 
hadron phase transition fel have indicated that, for a 



range of physical parameters controlling the transition, 
these effects are present. It is thus natural to consider if 
similar effects are present during heavy-ion collisions || . 

Whenever pre-transitional phenomena are relevant, 
one should expect modifications from the usual homoge- 
neous nucleation scenario, which is based on the assump- 
tion that critical bubbles of the hadronic phase appear 
within a homogeneous background of the QGP phase. 
The dynamics of weakly first-order transitions will be 
sensitive to the amount of phase mixing at T c : for large 
phase mixing, above the so-called percolation threshold, 
the transition may proceed through percolation of the 
hadronic phase, while for small amounts of phase mix- 
ing, by the nucleation of critical bubbles in the (inhomo- 
geneous) background of isolated hadronic domains, which 
grow as T drops below T c . An ideal quark gluon plasma 
in one dimension expands according to the Bjorken scal- 
ing, where T 3 t is constant fi"o|| . Assuming the initial tem- 
perature of the plasma produced at RHIC and LHC en- 
ergies to be 2 to 3 times T c , scaling implies that the time 
(At) taken by the plasma to cool from T\ to T c is of 
order a few fm/c, which could be comparable with the 
time scale of the subcritical hadronic fluctuations. On the 
other hand, the expansion rate of the early universe in 
the range T x < T < T c is slow enough |0Jl|] (At could 
be of the order of a few fi sees), that nonperturbative 
thermal fluctuations may achieve equilibrium. Another 
difference is that collisions at RHIC and LHC energies 
will lead to the formation of a highly (chemically) unsat- 
urated plasma, i.e., the initial gluon and quark contents 
of t he p lasma remain much below their equilibrium val- 
ues |l3|-|l5|] . A chemically unsaturated plasma will cool at 
an even faster rate than what is predicted from Bjorken 
scaling |ff^E8|. The cooling rate will also be accelerated 
further if expansion in three dimensions is considered. 
Therefore, we will show that although the equilibrium 
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density distribution of subcritical hadron bubbles is sig- 
nificant - particularly when the transition is weak - un- 
like the situation in cosmology, they do not contribute 
strongly to phase mixing. For the range of parameters 
we investigated, of relevance for RHIC and LHC energies, 
the plasma cools so rapidly that the subcritical bubbles 
do not have time to reach their equilibrium distribution 
and promote substantial phase mixing: significant super- 
cooling is possible in a (nearly) homogeneous quark gluon 
state. 



II. SUB-CRITICAL BUBBLE FORMALISM 

To study the dynamics of a first order phase transition, 
we use a generic form of the potential in terms of a real 
scalar order parameter <fi given by jl|J§, 



V{4>) = a(T)0 2 -bT(/) 3 + c4 4 



(1) 



The parameters a, b and c are determined from physical 
quantities, such as the surface tension (a) and the cor- 
relation length (£) of the fluctuations, and also from the 
requirement that the second minimum of the above po- 
tential should be equal to the pressure difference between 
the two phases || . The bag equation of state is used to 
calculate the pressure in the two phases. The potential 
V{4>) has a minimum at (j) — and a metastable second 
minimum at 



36T + y/9b 2 T 2 - 32ac 
8c 



(2) 



below T < T\, In the thin wall approximation [19], b, c 
and Ti can be written as |§, 
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where B is the bag constant and 



Vb((f> r , 
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16f(T c ) 
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(4) 



is the height of the degenerate barrier at T — T c or at 
a{T c ) — b 2 T 2 /Ac. A wide spectrum of first-order phase 
transitions, ranging from very weak to strong, can be 
studied by either changing a or £ or both. For example, 
for a fixed value of £, the strength of the transition is 
controlled by c, becoming very weak first order or second 
order when a — > 0. 

We follow Ref. || to obtain the equilibrium number 
density of subcritical bubbles. Let n(R, t) be the number 
density of bubbles with a radius between R and R + dR 
at time t that satisfies the Boltzmann equation 



fh_ M |n + (1 _ 7)ro _ 7r+ . 



(5) 



The first term on the right-hand side is the shrinking term 
with velocity v = dR/ dt. The term Tq is the rate per unit 
volume for the thermal nucleation of a bubble of radius R 
of phase <fr = (f> + (hadron phase) within the phase = 
(QGP phase). Similarly, T + is the corresponding rate of 
the phase = within the phase </> = </>+. The factor 
7 is defined as the volume fraction in the hadron phase. 
Assuming r ~ T + (= T) for a degenerate potential at 
T = T c , we write for the rate 



T = AT exp 



F(<t>+) 
T 



(G) 



where A is a constant of order unity. Using the Gaussian 
ansatz for subcritical configurations; 
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the free energy functional 



F(<f>) = 4tt / r 2 dr 
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can be written as Q 

F(<t>+) = a R + (3 R 3 , 

where 
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and 
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The equilibrium number density (no) of subcritical bub- 
bles is found by solving Eq. (||) with dn/dt = and im- 
posing the physical boundary condition n(R — > oo) = 0. 
Using 70 ~ 47ri? 3 no/3, we get a coupled equation for 70, 
which can be solved to get 



7o = 



1 + 21 



where 
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i? 3 r(i? ,(j>+)dR 



(12) 



(13) 



We will consider the statistically dominant fluctuations 
with R sa £ and estimate 70 integrating Eq. (13) from £ 
to 00. Neglecting the shrinking term in Eq. mJTthe time 
dependent solution of n(£,i) can be written as || 

nte,f)=«o(0[l-exp{-g(0*}]. (14) 

where q(£) = [(87r£ 3 /3)r] and n (£) = r(f)/g(£). Alter- 
natively, in term of 7, the above solution has the form 
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7fe*) = 7o(0[l-exp(-g i)], 



(15) 



where qo = (47r£ 3 /3)r/7o. The relaxation time r = Qq 1 
depends on two factors 70 and T out of which only the 70 
is affected by shrinking (if included). Since we know the 
complete solution of 70 that includes shrinking [Eq. (O)] , 
Eq. ( |l5| ) can also be used to estimate its time dependence. 
Note that the presence of a shrinking term in Eq. (|5|) 
results in a reduction of 70 and also in a faster relaxation 
process. 



III. RESULT AND DISCUSSION 

First we consider the slow evolution of the medium as 
in the case of early universe jy],[^) so that the equilib- 
rium scenario is applicable. Figure 1 shows the plot of 
70 as a function of a at T = T c for a few typical val- 
ues of the prefactor A. We have fixed £ at 0.5 fm, T c 
at 160 MeV and v — l/v3. As expected, the equilib- 
rium hadronic fraction increases with decreasing a and 
becomes as large as 0.5 for weak transitions. Recent lat- 
tice QCD predictions suggest that the quark-hadron 
phase transition could be weakly first order with er values 
in the range 2-10 MeV/fm 2 . Therefore, the choice of a 
in the above range and A *~ 1 ]l9| would imply significant 
amount of phase mixing at T — T c so that homogeneous 
nucleation becomes inapplicable M . 
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can estimate the time At taken by the plasma to cool 
from T\ to T c as 
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FIG. 1. 70 versus a at T — Tc for a few typical values of 
A. £ is fixed at 0.5 fm and T c at 160 MeV. 

Next we consider the plasma expected to be formed 
at RHIC and LHC energies. Since the expansion of such 
plasma is much faster compared to the plasma at the 
early universe, it will be interesting to know the amount 
of phase mixing (the value of 7) built up by the time the 
plasma cools from T\ to T c . Assuming ideal scaling, we 
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where v = 3 in (1+1) dimensions. Since T\ depends on cr 
[see Eq. (||)], At will also depend on cr, being smaller the 
weaker the transition. In the standard scenario, we can 
assume the initial temperature To « 320 MeV and the 
formation time to ~ 1 fm. However, several perturbative- 
inspired QCD models Jl4|-^6[ suggest a very different col- 
lision scenario at RHIC and LHC energies, which lead to 
the formation of unsaturated plasma with high gluon con- 
tent. Such a plasma will attain thermal equilibrium in 
a short time to ~ 0.3 — 0.7 fm, but will remain far from 
chemical equilibrium. Since the initial plasma is gluon 
rich, more quark and anti-quark pairs will be needed in 
order to achieve chemical equilibration. The dynami- 
cal evolution of the plasma undergoing chemical equi- 
libration was studied initially by Biro et. al. Jl7| and 
subsequently by many others [[18[ by solving the hydro- 
dynamical equations along with a set of rate equations 
governing chemical equilibration. It was found that a 
chemically unsaturated plasma cools faster than what is 
predicted by Bjorken scaling, since additional energy is 
consumed during chemical equilibrium. Following Ref. 
]l8| , we have studied chemical equilibration and dynam- 
ical evolution of the QGP with two sets of initial con- 
ditions, HIJING H] and Self Screened Parton Cascade 
Model (SSPM) fjj^ as listed in table I. The Perturba- 
tive QCD inspired models like Parton Cascade Model 
(PCM) pi and HIJING (Heavy Ion Jet Interaction Gen- 
erator) WM arc generally used to simulate the nuclear 
collisions at collider energies on the level of microscopic 
parton dynamics. The PCM calculations describe the 
space time evolution of quark and gluon distributions by 
Monte Carlo simulations of relativistic transport equa- 
tions. The HIJING model also incorporates the pertur- 
bative QCD approach and multiple minijet productions, 
however, it does not incorporate a direct space time de- 
scription. Early PCM calculations were done by assum- 
ing a pt cutoff to ensure the applicability of the pertur- 
bative expansion of the QCD scattering process. In the 
recently formulated self screened parton cascade model 
(SSPM) [3], early hard scattering produces a medium 
which screens the longer range color fields associated with 
softer interactions. The screening occurs on a length scale 
where perturbative QCD still applies and the divergent 
cross sections in the calculation of the parton produc- 
tion can be regulated self-consistently without an ad hoc 
cutoff parameter. The numerical studies based on the 
parton cascade model suggest that the parton plasma 
produced in the central region is essentially a hot gluon 
plasma and the dynamics is mostly dominated by glu- 
ons. Gluons thermalize rapidly reaching approximately 
isotropic momentum distributions in a very short time 
scale. The densities of quarks and antiquarks stay well 
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below the gluon density and can not build up to the full 
equilibration values required for an ideal chemical mix- 
ture of gluons and quarks. The similar conclusions have 
also been drawn from the calculations based on HIJING 
approach. Though both PCM and HIJING are QCD 
inspired models, the two still differ in quantitative pre- 
dictions possibly due to different treatment of multiple 
parton interactions and collective effects. In the follow- 
ing, we take the initial conditions obtained both from 
HIJING and SSPM calculations at the time when parton 
momentum distribution becomes isotropic. We consider 
two dominant reaction channels qq ^ gg and gg ^ ggg 
that contribute to the chemical equilibrium. The fugac- 
ity \ g (q)(< 1) gives the measure of the deviation of the 
gluon (quark) density from the equilibrium value; chem- 
ical equilibrium is achieved when Aj's — > 1. For a detail 
discussion on chemical equilibration, we further refer to 



TABLE I. Initial conditions are taken from Ref. [20] as 
predicted by SSPM and HIJING calculations. The fugacities 
Ai's give a measure of the deviation of the gluon or quark 
densities from the equilibrium values. 



CODE 


ENERGY 


U (fm/c) 


Ti (GeV) 


X g 


A, 


V 


SSPM 


RHIC 


0.25 


0.668 


0.34 


0.064 


2.2 


SSPM 


LHC 


0.25 


1.02 


0.43 


0.082 


2.2 


HIJING 


RHIC 


0.7 


0.55 


0.05 


0.008 


1.9 


HIJING 


LHC 


0.5 


0.82 


0.124 


0.02 


1.8 



Figure 2 shows a typical example of the effect of chemi- 
cal equilibration on the cooling rate for SSPM initial con- 
ditions at RHIC energy (A g o = 0.34, A g o = 0.064, to = 
0.25 fm and T = 0.668 GeV). The dotted curve (marked 
as Tb) shows the cooling rate as a function of time which 
obeys Bjorken's scaling (T 3 t=constant) corresponding to 
the case of an equilibrated plasma (A s = X q = 1.0). In 
case of a chemically unsaturated plasma for which the 
values of initial fugacities are much less than unity, the 
hydrodynamical expansion of the plasma proceeds along 
with chemical equilibration. As a result, both X g and 
X q increase with time as well as the temperature (shown 
by dashed curve) drops at a faster rate as compared to 
the Bjorken's scaling. The solid circles show the tem- 
perature given by (T v inconstant for v = 2.2). In this 
work, since we are interested only in the cooling rate, we 
skip the details of the calculation and parameterize the 
cooling rate in terms of v in the range T\ < T < T c (i.e. 
T^f^const). In table I, v has been listed for two sets 
of initial conditions obtained using HIJING and SSPM 
models at RHIC and LHC energies. Note that v < 3 
implies a faster cooling. Figure 3 shows the plot of At 
as a function of a as obtained from Eq. (16) for differ- 
ent v values. The time At depends on the initial values 
of the temperature To, formation time to and also on the 
cooling rate v. However, except for the SSPM initial con- 
ditions at LHC energies, values of At obtained with other 
initial conditions have nearly similar values. 
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FIG. 2. The temperature T and fugacity A as a function 
of time t. The description of the various curves are given in 
the text. 
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FIG. 3. At as a function of a for various initial conditions 
as shown in table I. 



Next we proceed to estimate the density of subcriti- 
cal hadron bubbles built up at t = At. Figure 4 shows 
7W/7o as a function oft at three different a values. The 
equilibration rate of the subcritical hadron bubbles of a 
given radius depends on the ratio r/70. Although both 
r and 70 are larger for weaker transitions, their ratio 
decreases with decreasing a. Therefore, as can be seen, 
equilibration is faster for a stronger transition as com- 
pared to the weak one. 
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FIG. 4. The ratio 7(f) /70 as a function oft at three typical 
values of a for j4=1. 

Figure 5 shows the fraction of the density built up at 
time t — At as a function of a with different initial con- 



ditions. Although the equilibrium density distribution of 
subcritical hadron bubbles increases with decreasing a, 
the time At decreases with decreasing a. As a result of 
these two competing effects, 7 at t = At shows a peak 
at around a m 20 MeV/fm 2 . The equilibrium fraction 70 
depends on the ratio A/v, which increases either due to 
increase in A or decrease in v. However, the variation in 
A and v act differently on qo as the nucleation rate F de- 
pends only on A. Therefore, we study the effect of A and 
V on 70 and 7 separately. Figure 6(a) shows the plot of 
70 (upper curves) and j(t) (lower curves) as a function of 
a at ^4=5, 10 and 20 respectively. Other parameters are 
v = 1/V3, T c = 160 MeV and £=0.5 f m . As expected, 
70 goes up as A increases. The increase in 70 for A from 
5 to 20 is about 1.5 to 2 times, but the nucleation rate T 
goes up by a factor of 4. Therefore, the ratio r/70 also 
goes up resulting in a faster equilibrium. The net conse- 
quence is both 70 and j(t) go up with increasing A. For 
the calculation of j(t), we have used SSPM and RHIC 
initial conditions. Further, we would like to mention here 
that although we have varied A up to 20, the value of A 
more than unity is unrealistic. A recent work by us p3] 
and also studies in ref p3] suggest A « 1. However, 
the ratio A/v can also go up with decrease in v, which 
we study in figure 6(b). Figure 6(b) shows 70 and 7 for 
v = c = 1 (upper limit), 1/4 and 1/12. This corresponds 
to a A/v ratio of 5, 20 and 60 respectively. Therefore, 70 
goes up with decreasing v as expected. Since A is fixed, 
r does not change, but q$ decreases with increasing 70 
resulting in slower equilibration. As a result, 7(f) does 
not build up at all. It is also interesting to note that j(t) 
does not get affected much by the choice of v although 
70 has a strong dependence on it. The j(t) only depends 
on parameter A. This aspect is interesting. 
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FIG. 5. The fraction 7 at t = At as a function of a at 
A = l and v = 0.577. 
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FIG. 6. (a) The fraction 70 (upper curve) and 7 at t = At 
(lower curve) as a function of a at A — 5, 10, 20 and v — 0.577. 
(b) Same as above at A— 5, but for different v values with 
SSPM initial conditions at RHIC energy. 

From the above studies (from figures 5 and 6), we can 
conclude that, the fraction in the range 2MeV/fm < 
a < lOMeV/fm does not build up to a significant level 
due to rapid cooling of the plasma, although the equilib- 
rium concentration is fairly large. It may be mentioned 
here that we have considered expansion only in (1+1) di- 
mensions. Inclusion of transverse expansion, significant 
at RHIC and LHC energies, will accelerate the cooling 
rate further, reducing the amount of phase mixing con- 
siderably. Since phase mixing at T = T c is negligible, 
the plasma will supercool and the phase transition may 
proceed by the nucleation of critical-size hadron bub- 
bles within a (nearly) homogeneous background of the 
metastable QGP phase. 

We have also studied the effect of other parameters like 
T c and £ on 7. Figure 7(a) shows the plots for various T c 
values at A=5. The nucleation rate decreases with de- 
creasing Tc [see Eq. (||)] resulting in a decrease in 70. On 
the other hand, smaller T c will result in larger At, which 
may increase "f(t). However, as shown in figure 7(a), the 
variation in j(t) with T c is not very significant although 
70 depends on it. Similarly, figure 7(b) shows the plots 
at various £ (0.5, 1.0 and 2.0). Increasing £ suppresses 
70 and j(t) particularly when the transition is strong. 
Therefore, the effect of other parameters like v, T c and 
£ on 7 are not very significant to promote phase mixing. 
The prefactor A is the only sensitive parameter on which 
j (At) depends. While the choice of A ss 1 is quite rea- 



sonable JlSt] , we have also varied A from 1 to 20 and did 
not find significant phase mixing particularly when a is 
small. 
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FIG. 7. (a) The fraction 70 (upper curve) and 7 at t = At 
(lower curve) as a function of a for different values of Tc 
at A — 5 and v = 0.577. (b) Same as above at A— 5, but 
for different £ values with SSPM initial conditions at RHIC 
energy. 



IV. CONCLUSION 

In conclusion, we have studied the effect of phase mix- 
ing promoted by thermal subcritical hadron bubbles dur- 
ing a first-order quark-hadron phase transition as pre- 
dicted to occur during heavy-ion collisions. Although the 
equilibrium density distribution of these subcritical bub- 
bles can be quite large, their equilibration time-scale is 
larger than the cooling time-scale for the QGP. As a con- 
sequence, for RHIC and LHC energies, they will not build 
up to a level capable of modifying the predictions from 
homogeneous nucleation theory. The phase transition 
may proceed either through the nucleation of critical size 
hadron bubbles in a (nearly) homogeneous background of 
the supercooled quark-gluon plasma or through spinodal 
decomposition if nucleation rate is not significant 
This situation is to be contrasted with the cosmological 
quark-hadron transition, where substantial phase mixing 
may occur, altering the dynamics of the phase transi- 
tion. We would also like to add here that even though 
our calculations rule out the role of subcritical bubbles, it 
is possible that impurities may increase the decay time- 
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scale and no real supercooling will be measured, as is the 
case with many condensed matter systems. The ques- 
tion, however, remains as to what these impurities, if any, 
might be in this context. One possibility - ruled in this 
work - is that the subcritical bubbles, being seeds for 
nucleation, may act as impurities [^4| . However, other 
possibilities, as the presence of condensates, may exist 
and should be considered in the near future. If there is 
supercooling there will be an extra entropy production 
which will reflect on the final hadron multiplicities. In 
this case, subcritical bubbles are not present, or are ir- 
relevant. On the other hand, if the transition is first order 
and no extra entropy is observed, subcritical bubbles (or 
unknown impurities...) do play a role. 
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